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Abstract 

We report on single-cluster Monte Carlo simulations of the Ising, 4-state 
Potts and 10-state Potts models on quenched ensembles of planar, tri-valent 
($ 3 ) random graphs. We confirm that the first-order phase transition of the 
10-state Potts model on regular 2D lattices is softened by the quenched con- 
nectivity disorder represented by the random graphs and that the exponents 
of the Ising and 4-state Potts models are altered from their regular lattice 
counterparts. The behaviour of spin models on such graphs is thus more 
analogous to models with quenched bond disorder than to Poisonnian ran- 
dom lattices, where regular lattice critical behaviour persists. 

Using a wide variety of estimators we measure the critical exponents for 
all three models, and compare the exponents with predictions derived from 
taking a quenched limit in the KPZ formula for the Ising and 4-state Potts 
models. Earlier simulations suggested that the measured values for the 10- 
state Potts model were very close to the predicted quenched exponents of the 
four-state Potts models. The analysis here, which employs a much greater 
range of estimators and also benefits from greatly improved statistics, still 
supports these numerical values. 



1 Introduction 



There has recently been some interest, and no little controversy, regarding the critical behaviour 
of systems with quenched bond disorder in 2D Q. It has been known for some time that the first- 
order phase transition displayed by q > 4-state Potts models on regular lattices is softened by the 
introduction of the quenched bond disorder to a continuous transition ||] , though the universality 
class of this transition and its dependence on the strength and nature of the bond disorder are 
still not completely clear ||, ||]. Models which already display a continuous transition in the 
pure case appear to have their critical exponents altered by the bond disorder H ||, [|] provided 
that the critical exponent a purc of the specific heat is positive. In addition, a qualitatively new 
phenomenon in the form of multi-fractal scaling of local correlators has also been predicted [|| 
and measured ||] when quenched bond disorder is present. It should be remarked that all of 
the theoretical results in the bond disordered case are perturbative in nature, which is in large 
measure the source of the controversy surrounding the various predictions for critical exponents, 
since the domain of validity is unclear. 

Another type of disorder that might be imposed is quenched connectivity disorder. Very high 
accuracy numerical simulations have shown that spin models on Poisonnian random lattices in 
both 2D III]] and 3D [11] stay stubbornly identical to their regular lattice brethren - there is no 



sign of the effects observed with quenched bond disorder. However, a different picture emerges 
when one considers spin models living on a quenched ensemble of tri-valent ($ 3 ) planar graphs, 
as generated by simulations of 2D quantum gravity. In this case the connectivity disorder is 



sufficiently strong for q > 4 Potts model transitions to be softened to continuous transitions [12 



and for q < 4 exponents to be modified from their regular 2D lattice values [13]. In this respect 
such planar random graphs appear to be much more akin to random bond disorder models than 
to the Poisonnian random lattices considered in [|l^, [TT[] . One very interesting feature of the 
3> 3 graphs is that exact, rather than perturbative, predictions for exponents exist in the q < 4 
case, by virtue of taking a quenched limit H, 14] in the KPZ formula which gives the weights of 



conformal operators when they are coupled to 2D quantum gravity. 

In the quantum gravity and string theory context one is typically interested in annealed 
rather than quenched connectivity disorder, in which the lattices and spins are interacting on 
the same time scale, providing a discrete analogue of the back reaction in continuum theories of 
gravity. In this case the relation between the bare (A) and dressed (A) weights is given by the 
KPZ relation || 



A _ V1-C + 24A- 

In order to calculate the dressed weights in the quenched case one should take c = in the KPZ 
relation to get 

A - v/T+24A-l (o) 

'-^quenched , j \") 

which now gives non-rational weights. Indeed, Cardy || has recently pointed out that the nth 
power of a correlator with bare weight A averaged over the disorder will scale not as nA, but 
rather 



A„= ^™-' , (3) 
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with a "typical" Q value being governed by | 

ox 



dn 



= 3A. (4) 

n=0 



In this paper we concern ourselves exclusively with measuring quantities which correspond 
to n = 1 in the formula above, such as the specific heat, susceptibility and magnetisation. The 
picture suggested by the earlier q = 2 and q = 10 simulations in [|l^, |l3| was that measurements 
were in accordance with the predicted quenched exponents for q < 4 and that the exponents of 
the quenched q > 4 models were "stuck" at the q = 4 values, as evinced by the good agreement 
between the measured q = 10 exponents and the quenched q = 4 predictions. This was somewhat 
similar to the original scenario postulated for the quenched bond disordered Potts models, where 
it was suggested that all bond disordered q > 4 Potts models displayed Ising-like criticality. An 
obvious test of the scenario for quenched connectivity disorder is to perform simulations at 
q = 4, as well as improving the quality of the measurements for q = 2 and q = 10 in order to 
get sharper estimates of the exponents in these cases. This is precisely what we do here. 

In what follows we briefly describe the simulation methods, an extension of those used 
in [p^[ , before going on to discuss the analysis of the results and our best estimates for the 
various exponents for the Ising, 4-state Potts and 10-state Potts models. We finish with some 
observations on the values we obtain. 



2 Simulation and Measurements 

We used the standard definition of the g-state Potts model partition function and energy in all 
the simulations, 

^Potts = 5>-^; E = -Y J S* i *v °* = !>•••> 1, (5) 

Wi} {ij) 

where (5 = J/ksT is the inverse temperature in natural units, 5 is the Kronecker symbol, and 
(ij) denotes the nearest-neighbour bonds of random <I> 3 graphs (without tadpoles or self-energy 
bubbles) with N sites. In this study we considered the cases q = 2 and 4 (with N = 500, 
1000, 2 000, 3 000, 4 000, 5 000, and 10 000) which in the pure model exhibit second-order phase 
transitions, and the case q = 10 (with N = 250, 500, 1 000, 2 000, 3 000, 5 000, and 10 000) which 
in the pure model undergoes a first-order phase transition. 



The simulations were carried out using the Wolff single-cluster update algorithm [16]. For 



each lattice size we generated 64 independent graphs using the Tutte algorithm [17], and per- 
formed 500K equilibration sweeps followed by up to 10 million measurement sweeps in order to 
obtain 500K independent measurement sweeps for each lattice size. The runs were carried out at 
several f5 values near the transition point and time series of the energy E and the magnetisation^] 
M = (gmaxjnj} — N)/(q — 1) recorded for each graph. In what follows the per-site quantities 
are denoted by e = E/N and m = M/N, the thermal averages on each individual graph by 
(. . . ) and the quenched average over the different graphs by [. . . ] av . 

From the time series of e and m it is straightforward to compute in the finite-size scaling (FSS) 



region various quantities at nearby values of (3 by standard reweighting |18[ techniques. Some 



1 Rather than average. 

2 Where n; < N denotes the number of spins of "orientation" i = 1,.. .,q in one lattice configuration. 
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care must be taken with the reweighting range in the presence of quenched averaging, but we 
confirmed that direct measurements of both the susceptibility and specific heat from fluctuations 
and numerical derivatives were in accordance with the values deduced from reweighting in several 
representative cases. Comparisons for q = 10 and N = 2000 around the reweighting point of 
Po = 2.22 are shown in Figs. 1 and 2. 

To estimate the statistical (thermal) errors for each of the 64 realizations, the time-series data 
was split into bins, which were jack-knifed IS] to decrease the bias in the analysis of reweighted 
data. The final values are averages over the 64 realizations which will be denoted by square 
brackets [. . . ] av , and the error bars are computed from the fluctuations among the realizations. 
Note that these errors contain both the average thermal error for a given realization and the 
theoretical variance for infinitely accurate thermal averages which is caused by the variation 
over the random graphs. 

From the time series of the energy measurements we compute by reweighting the average 
energy, specific heat, and energetic fourth-order cumulant, 



u(/3) 
V(f3) 



[(E)U/N, 

P 1 N[{e 2 ) - (e) 2 ] 

[1 



(6) 
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Similarly, we derive from the magnetisation measurements the average magnetisation, suscepti- 
bility, and magnetic cumulants 



x(P) 
U 2 (P) 

U 4 ((3) 



[<H>]av, 

PN[(m 2 )-(\m\)\ v , 
(m 2 ) 



[1 
[1 



3(|m|) 2jav ' 
(m 4 ) 



(7) 



3(rn 2 ) 2 v 

Further useful mixed quantities involving both the energy and magnetisation are defined by 



d[{\m\)] t 
dp 

<*ln[(|m|)] t 
dp 

dln[(m 2 )] £ 
dp 



[(\m\E)-(\m\)(E)] ay , 
\m\E) 



(H) 

(m 2 E) 
(m 2 ) 



(E) 
(E) 



(8) 



The dynamical aspects of the simulations are encoded in the autocorrelation functions and 
the associated integrated autocorrelation times f. It is customary [ 20 1 to convert the f thus ob- 
tained by multiplying with a factor / = ngi p (\C\) /N to a scale where, on the average, measure- 
ments are taken after every spin has been flipped once (similar to, e.g., Metropolis simulations). 
For quenched, random systems this procedure is not unique due to the necessary average over 
realizations, since one can take either [r] av = [/ • f] av or [/] av • [f] av . The differences between 
the two averaging prescriptions turn out, however, to be extremely small in practice. 



See below for some subtleties in the ordering of averages in the cumulants. 
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One finds that the autocorrelation times for q = 2 stay roughly constant with increasing 
system size, being [r e ] av ~ 3 — 4 for the energy and [r m ] av ~ 1.6 — 2.2 for the magnetisation. For 
q = 4 scaling behaviour is visible with [r e ] av ~ 12 — 18 and [r m ] av ~ 7 — 10, giving a dynamical 
exponent z/D ~ 0.064(10) for the energy. Power law scaling is much more pronounced for the 
q = 10 model with [r e ] av ~ 60 — 500 and [x m ] av ~ 40 — 350 and much larger dynamical exponents 
for both the energy and magnetisation, z/D ~ O(l). 

The self- averaging properties of the ensemble can be investigated by considering the proba- 
bility density for the r's, -P(r), rather than the average values, [r ejm ] av . One would expect the 
cumulative distribution F(r) = P(t')cIt' to tend to a step function for increasing system size 
in a self-averaging system. This is observed not to be the case for all the models simulated, 
giving clear evidence of non-self-averaging behaviour. These observations can be put on a more 
quantitative basis by looking at data collapse with the scaled variable Ar/[r] av [21|, where At 
is the standard deviation. 



3 Data Analysis and Results 

In the infinite- volume limit the various measured quantities exhibit singularities at the transition 
point. In finite systems the singularities are smeared out and scale in the critical region according 
to 

C = C rcg + N a / vD f c (x)[l + ...], 
X = N'^ D f x (x)[l + ...], 

= N~^ D f m (x)[l + ...], (9) 

= N^-^ D f m> (x)[l + ...}, 



[{\m\)\ 
d[(H)] 



dp 

dln[(|m| p )] 
d0 



N^ D f p (x)[l + ...], 



= N^ D f U2p (x)[l + ...], 

where C reg is a regular background term, u, a, (3, and 7, are the usual critical exponents, and 
the fi(x) are various FSS functions with 

x = ((3- (5 c )N l '" D (10) 

being the scaling variable. [1 + . . . ] indicates correction terms which become unimportant for 
sufficiently large system sizes N. We have expressed the scaling relations in terms of the total 
number of vertices N rather than the linear size L since the fractal dimension D of the graphs is 
a priori unknown. Numerical simulations and various analytic approaches suggest that D = 4 
[^2j ] for the ensemble of graphs we are considering, but we shall not need this explicitly for our 
analysis here. By rearranging equ.(|i~0|) one finally obtains the standard scaling relation for the 
peak-locations (pseudo-critical points) (3 C {N) on finite graphs, 

(3 c (N) = p c + aN- 1 /» D , (11) 

with a being a constant. 
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One further issue of principle remains. In the presence of quenched disorder there are three 
equally plausible ways of defining the various cumulants. For instance, if we take the fourth 
order magnetic cumulant as an example we could define not only the form used above 




(i) 



1 - 




J av 



(12) 



but also the variants 





(3) 



1 - 



1 - 



[(^ 4 )L 

3[(m 2 ) 2 ] av 

[<™ 4 >L 

3[<m 2 )]L' 



(13) 



and the correct choice is not immediately clear. We can hedge our bets in the scaling analysis 
by including all three of the variants without prejudice in order to check their consistency. In 
simulations with poor statistics per realisation (such as typically in spin-glass studies) usually 

(3) 

C/4 is taken since with that choice (systematic) bias effects are minimised. In our case, however, 
the statistics for each realisation is so large that there is no reason to favour one definition over 
the other on technical grounds. 

Without further ado, we start the analysis by estimating the value of 1/vD. With the 
wealth of available estimators we have various tactics available for the extraction of \/vD. 
One possibility would be to use the maxima of each of dU2/dP, dll^/dp (all variants for both 
cumulants), <iln[(| ?7i|)] av /(i/3, and c?ln[(?n 2 )] av / df3 as pseudo-critical points and then evaluate 
the scaling of each of these quantities at their own maxima to extract the exponent. Another 
would be to evaluate the scaling of each of these quantities at all of the available pseudo-critical 
points. A global estimate is then extracted by performing a direct or error weighted average. 
In both these cases we take a fairly conservative estimate for the errors by using the smallest 
contributing error bar 

We present the results from both approaches in Table 1. In obtaining these estimates we 
have dropped the smallest graph sizes in all cases and used for definiteness as our definition 
of the cumulants for the results presented in the table: and give values that are 

indistinguishable within the error bars. The results are quite stable to the deletion of the next 
smallest size, but the quality of the fits declines somewhat when this is done. In all of the 
listed fits the quality of fit Q was very good, the lowest being » 0.3, with most being as high 
as 0.8 — 0.9. For comparison we have included the prediction of the quenched KPZ formula 
(equ.(2)), the standard KPZ exponents and the regular 2D lattice exponents in the lower box. 
Since the q = 10 model has a first-order transition on a regular 2D lattice there is no direct 
prediction in this case. 



4 The largest contributing error bar would certainly be too pessimistic; this choice probably errs on the side of 



caution too. 
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Table 1: Fit results for the critical exponent 1/uD. 



quantity 


type 


q — z 


Q — 4 


q — ru 


ATT IAR 
CLU2/ up 


at maximum 


n "io( 1 ^ 


n A r >( r >\ 


u.oy^o ) 






32(1 ) 


D A?\C? \ 


D 61 (1) 




weighted av. 


0.31(1) 


0.44(2 


0.60(3) 


dU 4 /d(3 


at maximum 


0.32(2) 


0.39(2) 


0.62(3) 




average 


0.31(1) 


0.40(2) 


0.59(2) 




weighted av. 


0.31(1) 


0.40(2) 


0.59(2) 


dhx[(\m\)U/dl3 


at maximum 


0.36(1) 


0.42(1) 


0.56(1) 




average 


0.36(1) 


0.43(1) 


0.56(1) 




weighted av. 


0.36(1) 


0.43(1) 


0.56(1) 


dln[(m 2 )] av /d(3 


at maximum 


0.36(1) 


0.42(1) 


0.56(1) 




average 


0.36(1) 


0.43(1) 


0.57(1) 




weighted av. 


0.36(1) 


0.43(1) 


0.57(1) 


meta-average 




0.34(1) 


0.42(1) 


0.58(1) 


Quenched 




0.3486 . . . 


0.589 . . . 




KPZ 




0.3333 . . . 


0.5 




Regular 




0.5 


0.75 





Looking at the results in Table 1 it is clear that the estimates of 1/vD are not consistent with 
those for regular 2D lattices, giving a clear indication that the planar random graphs are different 
in this respect from Poisonnian random lattices. The average of averages or "meta-average" for 
q = 2 is compatible with both the quenched and KPZ values at the level of accuracy we have 
achieved, but that for q = 4 matches none of the possible theoretical predictions. Remarkably, 
the estimated q = 10 values are a good fit to the quenched q = 4 predictions, as we have already 
noted in [12|, and the numerous additional estimators here add extra weight to this observation. 
The difference between the theoretical quenched and KPZ values for q = 4 is sufficient for the 
q = 10 estimates not to be consistent with the theoretical q = 4 KPZ value. 

It is also noteworthy that the q = 10 measurements (and also the q = 4 quenched theory 
predictions) violate a supposedly general bound derived by Chayes et al. [23] for quenched 
systems, uD > 2, since vD ~ 1.72(3) from the q = 10 measurements. Hyperscaling, a/uD = 
2/vD — 1, implies that ct/vD should be negative if the bound holds, but we find below that 
direct fits to the specific heat for q = 10 also give a positive value (0.21(1)) that is compatible 
with that deduced from hyperscaling (0.16(1)). The measured values of vD for q = 2 and 
q = 4, on the other hand, are consistent with the bound. Whether the failure of the q = 10 
model to observe the bound is a consequence of the technical details of the averaging preocedure 
as suggested in [24] or a result of long range correlations in the disorder (which is due to the 
curvature correlations for the Liouville action in the original 2D gravity theory used to generate 
the graphs) is unclear. 

We now use our best estimates of 1/vD to extract the critical coupling f3 c by performing a 
linear two-parameter fit, equ.(ll), using the maxima of C and x along with derivatives of the 
three variants of U2 and U4 and derivatives and logarithmic derivatives of the magnetization 
as estimators of the pseudocritical points. This gives the eleven estimators which are shown in 
Table 2. A global estimate is again extracted by performing both a straightforward average and 
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an error weighted average. It is noteworthy that the estimated critical couplings are compatible 
with those found in simulations of the models coupled to 2D quantum gravity pi i.e., with 



annealed rather than quenched connectivity disorder. This was already remarked on in Jl^, 13] 
where non- linear three parameter fits were employed. 



Table 2: Fit results for the pseudocritical couplings (3 C 



at Pm&x °f 


q = 2 


q = 4 


q = 10 


c 


1.539(5) 


1.836(1) 


2.244(1) 


X 


1.562(6) 


1.834(2) 


2.246(1) 


dU^/dp 


1.551(5) 


1.831(1) 


2.244(1) 


dU^ 2) /dp 


1.550(5) 


1.834(1) 


2.242(1) 


dU^ 3) /dp 


1.550(5) 


1.831(2) 


2.243(1) 


dU^/dp 


1.558(7) 


1.832(2) 


2.244(1) 


dU^/dp 


1.569(4) 


1.841(2) 


2.240(2) 


duf ] /dp 


1.562(5) 


1.834(2) 


2.242(1) 


d(m)/dp 


1.561(3) 


1.834(1) 


2.245(1) 


d log (m)/dp 


1.562(5) 


1.837(2) 


2.243(1) 


dlog(m 2 )/dp 


1.570(5) 


1.838(1) 


2.243(1) 


average 


1.558(3) 


1.835(1) 


2.244(1) 


weighted av. 


1.558(3) 


1.835(1) 


2.244(1) 



The crossing points of the various definitions of the fourth-order cumulant in eqs. (@) and 
(HH) provide further estimates of p c . The error bars from these measurements are much larger 
than those in Table 2 due to the spread of crossings for different system sizes, but the estimates 



from l/[ , and U^' are all consistent, giving f3 c = 1.58(2) for q = 2, /3 C = 1.85(2) for 

q = 4, and P c = 2.244(4) for q = 10. The rather wide spread in the estimated p c is reflected by 

(2) 

a similar spread in the values of t/4 at P c , usually denoted by U*. The variant U\ appears to 
be the best behaved for all q giving U* ~ 0.55(3). 

Moving on to a/vD, the Fisher hyperscaling relation in the form 



r(3) 



a 



2 



1 



(14) 



gives a/vD = —0.32(1), —0.16(1), and 0.16(1) for q = 2, 4, and 10, respectively, from the 1/vD 
values in Table 1, whereas the quenched predictions are a/vD = —0.303. . . and 0.177. . . for 
q = 2 and 4. Direct non- linear fits to a/uD using the scaling form C = C reg + dN a / uD for 
the maxima of the specific heat are unstable for q = 2,4 but give a/uD = 0.22(7) for q = 10. 
The constant in this fit is consistent with zero, so we carried out a log-log fit at all of the 
pseudo-critical points as shown in Table 3, giving a final average of a/vD = 0.21(1) for q = 10. 
Two representative series of points for the specific heat evaluated at its own maximum and the 
maximum of the susceptibility on the different graph sizes are shown in Fig. 3. This independent 
measurement of a/vD when q = 10 is thus still in reasonable agreement with the theoretical 
quenched exponent values for q = 4. 



5 On corresponding ensembles of graphs without tadpoles and self-energy bubbles. 
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Table 3: Fit results for the critical exponent a/vD for q = 10. 



at /? max of 


q = lQ 


C 


0.21(1) 


X 


0.21(2) 


dU ( 2 > /dp 


0.22(1) 


dU^/dp 


0.21(1) 


dU {3) /dp 


0.22(1) 


du[ 1] /dp 


0.21(1) 


duf ] /dp 


0.18(2) 


dui 3) /dp 


0.21(1) 


d(m)/dp 


0.21(1) 


d log (in) /dp 


0.21(1) 


d\og(m 2 )/dp 


0.22(1) 


average 


0.21(1) 


weighted av. 


0.21(1) 


q = 4 Quenched 


0.177 ... 


q = 4 KPZ 





q = 4 Regular 


0.5 



In Tables 4, 5, and 6 we list the measured values of the magnetic exponents j/vD, P/vD 
and (1 — P)/vD for q = 2,4, and 10. We have adopted a similar approach to that used in the 
estimation of 1/vD. Data from all but the smallest graph size was included in the fits, and the 
stability of the fits checked against the deletion of the next smallest sizes, proving in all the 
listed cases to be reasonable. The appropriate scaling relations in equ. (||) are evaluated at all 
the pseudo-critical points reported in the tables and a final averaged value and error weighted 
average calculated. For the blank entries in the tables no stable fit proved possible. In the 
tables the number of terms used in the averages for each exponent are listed in the average 
row, and various theoretical predictions for the exponents are listed for comparison. Since it is 
already clear that the q = 10 measurements bear a considerable similarity to the quenched q = 4 
predictions we have repeated the q = 4 theoretical values in the q = 10 table. The exponent 
ratio j/uD is evaluated from the scaling of the susceptibility Xi P/vD from the scaling of the 
derivative of the magnetisation and (1 — P)/vD from the scaling of the logarithmic derivative. 

The scaling relation 

4^ = 1-2^ (15) 

relates the two exponents in these tables. The directly measured values are all in reasonable 
agreement with those derived by using either exponent as input in this scaling relation. There is 
an accidental equality between the theoretical quenched values of j/i/D (and hence via equ.(|i~5|) 
of P/uD) for the Ising and 4-state Potts models which is not, however, reflected in the estimates. 
There is a steady decrease in ^/uD and corresponding increase in P/vD as q is increased, so 
the q = 4 and q = 10 measurements are clearly different from those for q = 2. There is still a 
striking agreement between the quenched q = 4 predictions and the measurements at q = 10. It 
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is also clear that the q = 4 measurements are in definite disagreement with the quenched q = 4 
predictions. 

In Fig. 4 we plot the data points and fits for the susceptibility in the Ising model. The 
individual points for x evaluated at its own maximum, the maximum of the specific heat, the 
maxima of the derivatives of the second and fourth order cumulants and the maximum of the 
derivative of the magnetisation are shown explicitly []. The fits from the different series are all 
in good agreement, as indicated by the tables. We also plot \ evaluated at its own maximum 
for q = 2, 4, and 10 in Fig. 5 in order to show the general trend in the exponent. 

4 Conclusions 

The qualitative conclusions of our extensive analyses are quite clear: the quenched bond disorder 
of the $ 3 graphs does alter the exponents of models which already possess a continuous transition 
on a regular lattice, as well as softening the first-order transition of the q = 10 model to a 
continuous transition. A quenched ensemble of $ 3 graphs with connectivity disorder thus shares 
many of the features of a system with quenched bond disorder. As noted in the introduction, 
other ensembles with quenched connectivity disorder such as Poisonnian random lattices do not 
- being more akin to regular lattices. 

At a quantitative level, however, the current batch of simulations pose rather more questions 
than they answer: the working hypothesis of the veracity of the quenched exponents is at best 
only weakly supported by the results. For the Ising (q = 2) model, the estimated value of 1/uD, 
0.34(1), is consistent with both the quenched and KPZ predictions. Although the estimates 
for the magnetic exponents [j/uD, (3/uD, (1 - (3)/vD] ~ [0.79(1), 0.11(1), 0.26(1)] are closer 
to the quenched [0.7094, 0.1452, 0.2033] than the Onsager [0.875, 0.0625, 0.4375] or KPZ 
[0.666.., 0.166..., 0.166...] values, any agreement is less than convincing. In mitigation, it is fair 
to point out that we have struggled in the past to obtain good estimates of the magnetic KPZ 
exponents on dynamical <I> 3 graphs without tadpoles and self-energy bubbles |f25f| , which is the 
class of graph we have used in the simulations here. This is probably due to large corrections 
to scaling, since including degenerate graphs (self-energy bubbles and tadpoles) appears to give 
a faster approach to the continuum limit p6fl . 

Taken en masse the estimated exponents of the q = 4 model, [1/uD, "y/uD, (3/vD, (1 — 
0)/vD] ~ [0.42(1), 0.75(1), 0.12(1), 0.34(1)] fit neither the quenched [0.589 ... , 0.7094 
0.1452..., 0.4433 ...], KPZ [0.5, 0.5, 0.25, 0.25] nor regular lattice [0.75, 0.875, 0.0625, 
0.6875] predictions, although one could argue that j/vD ~ 0.75(1) on its own is actually closer 
to the quenched prediction than the q = 2 model. It is possible that q = 4, which is subject to 
logarithmic corrections in both the regular lattice and KPZ cases, may require similar treatment 
in the quenched case but without more input on the form of these corrections fitting would be 
a futile exercise. 

The q = 10 model, on the other hand, provides us with a set of estimated exponents 
[1/vD, y/vD, P/vD, (1 - 0)/vD] ~ [0.58(1), 0.71(1), 0.12(1), 0.43(2)] which match the 
predicted quenched q = 4 exponents [0.589.., 0.7094.., 0.1452.., 0.4433..] extremely well. These 
estimated q = 10 exponents (and theoretical q = 4 values) hence violate the bound vD > 2 of 
|p3f| , which merits an explanation in its own right. 

6 We have shown only the cumulant results to avoid cluttering the graph, the alternative definitions give 
effectively identical results. \ evaluated at the maxima of the derivatives of the log of the magnetisation and the 
modulus squared of the magnetisation has also been dropped for clarity. The graphs for q = 4 and q = 10 show 
similar features and are not reproduced here. 
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It is difficult to quibble with the thermal statistics from the very long time series in the current 
batch of simulations, but one might find fault with the relatively modest number of replicas used 
in the disorder averaging. The replica to replica variation of the measured quantities does not, 
however, appear to have significantly skewed the measured values and error estimates. Looked 
at without theoretical prejudice the measured exponents suggest a slow variation with q that is 
akin to that observed in the quenched bond disorder simulations of [||, rather than values which 
stick at q = 4 and change no further with increasing q. It is thus possible that the agreement 
of the measured q = 10 exponents with the predicted q = 4 values is accidental. A simple way 
to settle this issue would be to simulate other values of q to investigate the variation, if any, of 
the exponents, particularly for q > 4. 

There is a second parameter that one can vary in such simulations, namely the central 
charge used in generating the graphs of the quenched ensemble. Since the partition function Zjv 
obtained on integrating out d scalar fields on $ 3 graphs with N vertices is 

z N = ( det c G)~ d/2 > ( 16 ) 

Geg(A) 

where G(A) is the class of graph being summed over and Cq is the adjacency matrix of the 
graph G, 

Qi if i = j, 

C G = { 

—riij if i and j are adjacent, (17) 

otherwise, 

one can use equ.fllBD to generate an ensemble of graphs to which one can associate a central 
charge d. 

Quenched simulations may then be carried out on this ensemble rather than the d = c = 
graphs used here. The appropriate d can be substituted into the KPZ formula to obtain 



predictions for the exponents in this case. This brings one to another puzzle: in [27| very 
good agreement was found between the predicted exponents for the Ising model on a quenched 
ensemble of graphs with d = — 5 and measurements. In this light, it is surprising that the 
agreement here for the Ising model on much larger d = graphs with better statistics is poorer. 
One could speculate that the scaling behaviour of the models improved as they became more 
"classical", i.e. as the effects of gravity were switched off (d — ► — oo). 

Various other aspects of the behaviour of spin models on quenched random gravity graphs 
that have only been touched on here merit further investigation. The clear evidence of non- 
self- averaging behaviour for all q and the autocorrelation scaling techniques used to quantify it 



are described in more detail in a companion paper [21]. Similarly, the multifractal scaling of 
spin correlation functions can also be investigated. Finally, as we have already noted, further 
simulations for other q (and also d) values would help to determine whether the quenched 
exponents were correctly describing the behaviour of the models and cast further light on the 
remarkable (accidental?) agreement between q = 10 measurements and the predicted q = 4 
exponents. 

In summary, spin models on $ 3 random graphs offer a useful framework for the exploration 
of quenched disorder and may even offer some advantages over bond disordered models given the 
availability of various exact, rather than perturbative, predictions for exponents. The results 
described here suggest numerous avenues for future work. 
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Table 4: q=2 Potts (Ising) fit results for the critical exponents j/uD, fi/vD, and (1 — (3)/vD. 



at /J max oi 


i i ~i 
j/vD 


Q 


p/vD 


Q 


(l-p)/vD 


Q 




U. IO{L ) 


n 1 a 
U.14 


n i a { o\ 

U.14(Z ) 




\J.l6\l ) 


n a n 
U.4U 


X 


0.79(1) 


0.13 


0.08(1) 


0.21 


0.26(1) 


0.78 


rlTT^ 1 AR 
CLU 4 j dp 


0.79(2) 


0.66 


0.12(1) 


0.42 


0.26(2) 


0.33 


ail 4 /dp 














dJjf ] /d(3 


0.80(2) 


0.46 


0.10(1) 


0.14 


0.28(1) 


0.13 


dU$ ' /dp 


0.77(1) 


0.75 


0.12(1) 


0.82 


0.25(1) 


0.76 


duP/d/3 


0.77(1) 


0.73 


0.12(1) 


0.64 


0.25(1) 


0.63 


dU^/dp 


0.77(1) 


0.78 


0.12(1) 


0.84 


0.25(1) 


0.62 


d[(\m\)] av /dp 






0.09(1) 


0.38 


0.26(1) 


0.80 


dln[(\m\)U/dP 


0.80(1) 


0.55 


0.10(1) 


0.36 


0.28(1) 


0.37 


dln[(m 2 )U/dp 


0.83(2) 


0.15 






0.30(1) 


0.32 


average (9,9,10) 


0.79(1) 


0.11(1) 


0.26(1) 


weighted av. 


0.78(1) 


0.10(1) 


0.26(1) 


Quenched 


0.7094 . . . 


0.1452 ... 


0.2033 . . . 


KPZ 


0.6666 . . . 


0.1666 ... 


0.1666 ... 


Regular 


0.875 


0.0625 


0.4375 



Table 5: q=4 Potts fit results for the critical exponents j/i/D, (5/vD, and (1 — P)/uD. 



at P max of 



7/VD 



Q 



P/vD Q 



{l-P)/vD Q 



C 

X 

dU^/dp 
dU^ /dp 
dU^ /dp 
dU^/dp 
dU^ /dp 
duf ] /dp 
d[(\m\)] av /dp 
dln[(\m\)U/dp 
dln[(m 2 )U/dp 



0.75(1) 



0.13 



0.75(1) 
0.77(1) 
0.75(1) 
0.75(1) 



0.82 
0.42 

0.80 
0.44 



0.12(1) 

0.13(1) 
0.10(1) 



0.50 

0.26 
0.99 



0.34(2) 
0.33(1) 



0.33(1) 
0.34(1) 
0.33(1) 
0.34(1) 



0.61 
0.37 



0.70 
0.41 
0.66 
0.73 



average (5,3,6) 



0.75(1) 



0.12(1) 



0.34(1) 



weighted av. 



0.75(1) 



0.11(1) 



0.34(1) 



Quenched 



0.7094 



0.1452 



0.4433 



KPZ 



0.5 



0.25 



0.25 



Regular 



0.875 



0.0625 



0.6875 
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Table 6: q=10 Potts fit results for the critical exponents j/uD, P/vD, and (1 — (3)/vD. 



at /3 max of 


j/uD 


Q 


P/uD 


Q 


(1 - P)/uD 


Q 


C 


0.71(2) 


0.13 










X 


0.72(1) 


0.31 


0.10(1) 


0.39 


0.43(2) 


0.10 


dU^ /dp 


0.73(2) 


0.72 


0.11(1) 


0.23 


0.43(2) 


0.53 


du[ 2) /dp 














duf ] /dp 


0.69(2) 


0.30 






0.41(2) 


0.22 


c#7 2 (1) /dp 


0.72(2) 


0.61 


0.12(1) 


0.32 


0.43(2) 


0.19 


dU^/dp 










0.42(2) 


0.16 


dU^ /dp 


0.71(2) 


0.37 


0.14(1) 


0.29 


0.43(2) 


0.12 


d[(\m\)l v /dp 


0.72(2) 


0.25 


0.11(1) 


0.10 


0.43(2) 


0.11 


dln[(\m\)U/dp 


0.71(2) 


0.10 






0.43(2) 


0.10 


dln[(m 2 )l v /dp 










0.43(2) 


0.11 


average (8,5,9) 


0.71(1) 


0.12(1) 


0.43(2) 


weighted av. 


0.71(1) 


0.12(1) 


0.43(2) 


q = 4 Quenched 


0.7094 . . . 


0.1452 . . . 


0.4433 . . . 


q = A KPZ 


0.5 


0.25 


0.25 


q = 4 Regular 


0.875 


0.0625 


0.6875 
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q = 10, N = 2000, P =2.22 

64 repl., 500,000 + 12*500,000 



2.0 




reweighting 
fluctuations 
num. deriv. 



2.4 



2.5 



Figure 1: The specific heat calculated from reweighting about (3$ = 2.22, fluctuations and 
numerical differentiation of the energy for q = 10 and N = 2000 (which is typical). 
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q = 10, N=2000, P =2.22 

64 repl., 500,000 + 12*500,000 (without 0-factor) 
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Figure 2: The susceptibility calculated from reweighting about (3q = 2.22 and from fluctuations 
for q = 10 and N = 2000. 
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Specific Heat vs N for q=1 
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Figure 3: Two representative fits (from the eleven used) for the scaling of the specific heat 
evaluated at its own maximum and at the maximum of the susceptibility %. 
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Susceptibility vs N for q=2 
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Figure 4: Fits to x when q = 2 to obtain j/uD at the maxima of X , C, dU^/dp, du1 ] /d(3 
and dm/ d(5 are shown on a log- log scale. 
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Figure 5: Fits to x & t its own maximum for q = 2,4, and 10. 
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